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^ ■ 1. Introduction 



The Kerr spacetime is one of the most important exact solutions to the Einstein vacuum 
equations. Its relevance stems from the uniqueness theorems for black holes which state 
that under rather general conditions the Kerr spacetime is the only asymptotically flat, 
stationary, vacuum black hole solution — see e.g. the introduction of [19] for a critical 
review of the issue of black hole uniqueness and the involved assumptions. Although, the 
Kerr spacetime is very well understood from a spacetime perspective, the same cannot 
be said if one adopts a 3+1 point of view — which would be the case if one tries to 
numerically calculate the spacetime from some Cauchy initial data. 

As soon as one moves away from a 3+1 gauge which is adapted to the stationary 
and axial symmetries — which can occur in some applications, in particular in numerical 
ones — an analysis of the Kerr spacetime and initial data sets thereof becomes very 
complicated. The explicit nature of the Kerr solution makes it tempting to perform 
detailed calculations in order to, say, verify a particular property of the spacetime. This 
approach usually leads to very lengthy expressions which can be very hard to analyse. 
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In a sense one could say that exact solutions contain too much information. In this case 
it can be more convenient to adopt more abstract approaches and, accordingly, it may 
prove useful to have at hand a characterization of Kerr initial data. 

The question of providing an invariant characterization of initial data sets for the 
Schwarzschild spacetime has been addressed in [28l [16] . In particular, the analysis of [T6] 
provides an algorithmic characterization of Schwarzschild data. That is, a procedure is 
provided to verify whether a given initial data set for the Einstein field equations will 
render a development which is isometric to a portion of the Schwarzschild spacetime. 

One of the most important algebraic properties of the Kerr spacetime is that its 
Weyl tensor is of Petrov type D. The close relation between vacuum spacetimes with a 
Weyl tensor of Petrov type D and Killing spinors has been exploited in [17] to provide 
a characterization of initial data sets whose developments will be of Petrov type D. 
This characterization relies on one being able to decide whether a set of overdetermined 
partial differential equations has solutions for a given initial data set. Accordingly, such 
a characterization is not algorithmic. Although not explicitly stated in [T7J, from that 
analysis it should be possible to obtain a characterization of Kerr initial data by adding 
some global conditions. 

The characterization of initial data sets discussed in [28j EE E] nas followed the 
general strategy of starting from a given tensorial (respectively, spinorial) spacetime 
characterization of the class of spacetimes under consideration. Necessary conditions on 
the initial data set are obtained by performing a 3+1 decomposition of the spacetime 
characterization. Given a set of necessary conditions on the initial data, it is then 
natural to address the question of sufficiency This is, usually, the most challenging part 
of the procedure as one has to discuss the evolution of complicated tensorial objects. 
The idea behind this is to show that if the necessary conditions are satisfied on some 
open subset of the initial hypersurface, then one can — possibly, under some additional 
assumptions — recover the spacetime characterization on the development of the open 
subset on the initial hypersurface from which one started. 

In this article a particular characterization of Kerr initial data is addressed. Our 
starting point is a certain spacetime characterization of the Kerr spacetime given in 
[2T1 [22J. This characterization was developed with the aim of providing an alternative 
way of proving the uniqueness of Kerr spacetime among the class of stationary, 
asymptotically flat black holes. This expectation has been recently fulfilled in [19J, 
where a proof of the uniqueness of Kerr which does not make assumptions on the 
analyticity of the metric has been given. At the heart of the spacetime characterization 
given in [2TL [22] — cfr. theorem [3] — is a certain tensor, the Mars-Simon tensor, whose 
construction requires the existence of a timelike Killing vector. The Mars-Simon tensor 
is a spacetime version of the Simon tensor, a tensor, defined in the quotient manifold 
of a stationary spacetime, which characterizes the Kerr spacetime — see |27j . 

Following the general strategy for the construction of characterizations of initial 
data sets out of spacetime characterizations, necessary conditions for Kerr initial data 
are deduced from a 3+1 splitting of the Mars-Simon tensor. Accordingly, one assumes 
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that the spacetime one is working with has a timelike Killing vector. This requirement 
can be encoded in the initial data by requiring that the data has a suitable Killing 
initial data (KID) . The Mars-Simon tensor has the same symmetries as the Weyl tensor, 
and hence its 3+1 decomposition can be given in terms of its electric and magnetic 
parts. In order to discuss the propagation of the Mars-Simon tensor we make use 
of a framework for discussing the propagation tensorial fields using superenergy-type 
tensors — this framework has been discussed in e.g. [3]. It should be pointed out that 
the characterization discussed in this article is not algorithmic. That is, like the one 
for type D initial data discussed in [IT] it depends on being able to decide whether a 
certain overdetermined system of partial differential equations admits a solution. 

The article is structured as follows: in section [2] our main conventions are fixed and 
relevant aspects of the 3+1 formalism are discussed. Section [3] discusses the properties 
and causal propagation of Weyl candidates — i.e. tensors with the same symmetries of 
the Weyl tensor. Section H] is concerned with the properties of the Mars-Simon tensor. 
Section [5] discusses the causal propagation of the Mars-Simon tensor. Section [6] applies 
the previous discussion to the construction of a characterization of Kerr initial data. 
Our main result is provided in theorem [61 In section 16.11 we particularize to the case of 
the Schwarzschild spacetime, where, with the aim of the results of [16] it is possible to 
obtain an initial data characterization which is algorithmic. Some concluding remarks 
are provided in section [7J Finally, some technical details, too lengthy to be presented 
in the main text are presented in the appendix A. 

All the tensor computations of this paper have been performed with the software 
xAct [23] . xAct is a suite of MATHEMATICA packages which has among its many 
features the capability to efficiently canonicalize tensor expressions by the use of powerful 
algorithms based on permutation group theory [23] ■ Currently no other software, either 
free or commercial, is capable to handle the tensor computations needed in this paper. 

2. Preliminaries 

Let (M,g/j, u ) denote a smooth orientable spacetime. The following conventions will 
be used: plain Greek letters a, [3, 7, . . . denote abstract indices and boldface Latin 
characters a, b, c, . . . will be used for component indices. The signature of the metric 
tensor g^ v will be taken to be (— , +, +, +), while R^ U a P , R^u = R° ' mv and Wap^ denote, 
respectively the Riemann, Ricci and Weyl tensors of g^ u . The tensor r] a p au is the volume 
element which is used to define the Hodge dual of any antisymmetric tensor — denoted 
by attaching a star * to the tensor symbol. Sometimes we will need to work with 
complex tensors in which case the complex conjugation of a tensor is denoted by an 
overbar. The operator symbolizes the Lie derivative with respect to the vector field 
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2.1. The orthogonal splitting 

Let n M be a unit timelike vector, n^n^ = — 1 defined on A4. Then any tensor or tensorial 
expression can be decomposed with respect to n M and the way to achieve it is the essence 
of the orthogonal splitting (also known as 3+1 formalism) which is described in many 
places of the literature — see e.g. [HI [12]). We review the parts of this formalism needed 
in this work. The spatial metric is defined by h^ u = g^ u + n^n v and it has the algebraic 
properties /i^ = 3, h°h av = h^ u . We shall call a covariant tensor T ai _ am spatial with 
respect to h^ u if it is invariant under h^ u i.e. if 

h ai Pi ■ ■ ■ h am p m T av .. am = Tp v ..p m , (2.1) 

with the obvious generalization for any mixed tensor. This property implies that the 
inner contraction of n M with T ai ___ am (taken on any index) vanishes. The orthogonal 
splitting of a tensor expression consists in writing it as a sum of terms which are tensor 
products of the unit normal and spatial tensors of lesser degree — or the same degree in 
which case the unit normal is absent. 

In order to find the orthogonal splitting of expressions containing covariant 
derivatives we need to introduce the spatial derivative which is an operator whose 
action on any tensor field T 01 '""^ ^ , p, q G N is given by 

P Pl...f)q pi Pp Pl Pq p A (J\...(Jq \ ) 

From equation (I2.2p it is clear that D^T ai '" ap f3i ^ is spatial. 

The results just described hold for an arbitrary unit timelike vector n M but in our 
framework we only need to consider integrable timelike vectors which are characterized 
by the condition n^V^TV] = (Frobenius condition). In this case there exists a foliation 
of A4 such that the vector field is orthogonal to the leaves of the foliation. We shall 
denote by {St}, t G / C R, the family of leaves of this foliation and S is called the 
initial data hyper surface — it is assumed that G /. The tensor h^ u plays the role of the 
first fundamental form for any of the leaves while the symmetric tensor K^ u defined by 

can be identified with the second fundamental form. Combining the previous definition 
with the Frobenius condition we easily derive 

= -K„ v - n^A v , (2.4) 

where = n p V p n^ is the acceleration of n M . By using these quantities one is in 
principle able to work out the orthogonal splitting of any tensorial expression. We 
supply below the explicit results of calculations which will be needed repeatedly in the 
sequel. 

• Orthogonal splitting of the volume element: 

(2.5) 

Here e a/ g 7 = n M ?7 MQ/ 3 7 is the spatial volume element which is a fully antisymmetric 
spatial tensor. 
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• Orthogonal splitting of an antisymmetric tensor. If H pu is an arbitrary 
antisymmetric tensor then its orthogonal splitting takes the form 

H pv = 2P [v n p] - Q p e^, (2.6) 

where P p = H pu n u , Q p = H* u n u are clearly spatial tensors. It is straightforward 
from the previous equation to find an analogous formula for the orthogonal splitting 
of H*. If we work with the self-dual part of H pu , denoted by 7i pu and defined by 

H pv = H pv + i H* u , (2.7) 

then its orthogonal splitting takes the form 

T~Cfiv = i e pia ,V - 2V [ „n„ ] , V p = P p + i Q p . (2.8) 

• Orthogonal splitting of the Weyl tensor: 

W puXp = 2 (l p [ X E p]u - l u[x E p]p - n [x B p]r e T pu - n [p B u]r e T Xp ) , (2.9) 

where 

E Tp = W TUpX n v n\ B rp = W; vpX n»n\ (2.10) 

denote the Weyl tensor electric and magnetic parts respectively, and l pv = h pu + 
n p n u . The tensors E pu and B pv are symmetric, traceless, and spatial. Equation 
(12.91) also holds for any rank-4 tensor possessing the same algebraic properties as 
the Weyl tensor. Such tensors will be designated collectively as Weyl candidates 
—see section [3] We can also choose to work with the self-dual Weyl tensor W puXp 
which is given by 

w^x P = l(w, vXp + iw; vXp ), (2.ii) 

Using (12.91) and the analogous formula for the orthogonal splitting of W* x one 
can obtain an expression for the orthogonal splitting of W puXp 

yVpv P \ = 2{L[\£p]p + lp,\p£\\v + i SpXa^ifiS^ a + i e pua n[ p S x "), (2.12) 

where 

£pu = W ppuX n p n x = -{E pv - i B pu ), £ (pu) = S pu , £ p p = 0. (2.13) 

Again a similar formula to ( 12. 12ft holds for any Weyl candidate. 

• Let T p be any spatial vector. Then the orthogonal splitting of V p T u is 

V P T U = D P T U - n p {K pu T p + £ H T U ) - T p A p n p n u - n v K pp T p . (2.14) 
Note that since T p is spatial then L^T V is also spatial. 



2.2. The Cauchy problem 

We briefly review the standard formulation of the Cauchy problem for the vacuum 
Einstein equations. In this formulation one considers a 3-dimensional connected 
Riemannian manifold (£, hy) — we use small plain Latin letters i, j, k, . . . for the abstract 
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indices of tensors on this manifold — and an isometric embedding : £ — > A4 . The map 
(f) is an isometric embedding if dffidffi = hij where di^djcjf realises the pullback, 
<p*, of tensor fields from A4 to £. The metric h%j defines a unique affine connection Di 
(Levi-Civita connection) by means of the standard condition 

Djhik = 0. (2.15) 

The Riemann tensor of Di is r^jy and from it we define its Ricci tensor by r^ = r l u - 
and its scalar curvature r = r\ — in £ indices are raised and lowered with and its 
inverse h % K 

Theorem 1. Let (£, h^) be a Riemannian manifold and suppose that there exists a 
symmetric tensor field Kij on it which verifies the conditions (vacuum constraints) 

r • K- •- /\''/\, ( 0. (2.16a) 
LVKij - DiK = 0, (2.165) 

where K = K\ . Provided that and are suitably smooth there exists an isometric 
embedding <fi of £ into a globally hyperbolic, vacuum solution (Af, g^ v ) of the Einstein 
field equations. The set (£, hij, K^) is then called a vacuum initial data set and the 
spacetime (Af, g^ v ) is the data development. Furthermore the spacelike hypersurface 
</>(£) is a Cauchy hypersurface in M . 

A statement of this theorem containing precise regularity conditions on hij, K^ is 
formulated in reference [20J. 

Under the conditions of the theorem dj we may construct a foliation of M with n 11 
as the timelike unit vector which is orthogonal to the leaves and, in the above notation, 
set E = 0(E). It is then clear that di^djcjfh^ = h^. Other key properties are 

d^dj^K^ = K^, d^dj^ ■ ■ ■ d jq (j) Vq (DnT Ul _ Uq ) = Di(d Ul ^ ■ ■ ■ d Vq( p Vq T Ul .. Mq ).(2.17) 

Using these properties and the Ricci identity for n M one can show that 

Eij = d^dj^E^ = nj + KKij - K ik K h p (2.18a) 

Bij ee d^dj^B^ = e^ l D\ k K m . (2.18b) 

Another property which is needed later on is the following one: for any spatial tensor 
Pa!...a p , P £ N with respect to we have the property 

P ai ...a P \m = 0*(P*i...« P ) = 0- (2.19) 
3. Weyl candidates 

Fundamental for our discussion will be tensors of rank 4 having the same algebraic 
properties as the Weyl tensor. More precisely, one defines a Weyl candidate as any rank 
4 tensor C a p^ fulfilling the properties 

C[a/3]7<5 = Ca/3-yS = = 0, C a [/3^s] = 0. (3.1) 
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From flHUJ we easily deduce that *C al3lS = C*^ = (l/2)r7 7<5(TT C Q(g ,7T and indeed C*^ s 
is also a Weyl candidate. Given a Weyl candidate C a ^, its VKey/ current J af3l is defined 
by 

JpyS = VcC" 1 ^. (3-2) 

The tensor Jg 7< 5 is trace-free, antisymmetric in the last pair of indices and it has the 
property J\p>y8\ = 0. From (13.21) we deduce the following identities 

^\aPpi\ia> — ~^Va/3'raJ* CT fJ , u , V 'fiJ^jg = Rft ^ 'S]a pp,i (3-3) 

where J*° ' = Vu,v pX J cr p\/2- Equations (I3.2p - (I3.3I) have counterparts involving C* 
which we will not discuss. Also, if we combine (13.21) and (13.31) we get a wave equation 
for the Weyl candidate: 

Vr7a/i O D Q & r* np a Sri i O P a f~> I Of?" ^ C 1 

a V L,/3 7X p - —4K y [p Ly x ]8f3a - ^ x p ^p-yaS + ^ X }»Pl + Z±t /3 [p x]5 7 a 

+V XP x^J*\ + 2V [X J p]Pt (3.4) 

The previous considerations hold regardless of whether the Weyl candidate is real or 
complex. Indeed, from any real Weyl candidate C pvpa we may construct a complex Weyl 
candidate C pvpa by means of 

C^pa = -(C pupa + i C pupa ). (3-5) 



A Weyl candidate constructed in such a way is self- dual 

^pupa = 1 Cpvpa- (3-6) 

Reciprocally, if a complex Weyl candidate C pvpa fulfils ( 13. 6 j) then its real part Re(C pup \) 
and its imaginary part lm(C pup \) are also Weyl candidates related via 



Im(C^ pA ) = ~ VpXar3 Re(C«P), Re(C pupX ) = -~ VpXa0 Im(C«P). (3.7) 



3.1. The causal propagation of a Weyl candidate 

Besides the basic algebraic properties of a Weyl candidate just explained we need to 
introduce a further concept which will play a crucial role in the sequel. 

Definition 1. Let (Ai,g pu ) be a spacetime and consider a Weyl candidate C ai 3 pv defined 
in the whole of M. . We say that the Weyl candidate C a p pv propagates causally on M. 
if for any embedded spacelike hypersurface B C M. the condition C al 3 pu \s = implies 
C a p pu = onD(B). 

Given a set B, a point p G M. belongs to the future Cauchy development D + {B) if any 
past-inextendible causal curve containing p intersects B. There is a similar notion of 
D~{B) and the set D{B) = D+(B)UD-(B) is called the Cauchy development of B. It is 
well-known that the interior of D(B) is globally hyperbolic and since B is an embedded 

| Note that sometimes the opposite convention is followed in the literature whereby the complex 
conjugate of C pvp<J is given the name of "self-dual" . 



Kerr initial data 



8 



spacelike hypersurface then it will be a Cauchy hypersurface for D{B) — see e.g. [22] 
for an elementary introduction to these concepts of causality theory. 

The simplest example of a Weyl candidate which propagates causally is the Weyl 
tensor W a p^ v of a vacuum spacetime — see e.g. [5]. It is possible to address the causal 
propagation of a Weyl candidate by following the general techniques given in [4] — see 
also [26] — in which the causal propagation of any tensor is analysed. These ideas are 
reviewed next. 

An essential object to study the causal propagation of a Weyl candidate is its Bel- 
Robinson tensor. Given a real Weyl candidate C^s, we define its Bel-Robinson tensor 
by 

T a pj5 = C^^Cppy^ + CJ^^Cp^p — -g a pg-ysCp Up \C^ upX . (3.8) 

The previous definition has been described in a more general framework in [25] as the 
superenergy tensor of the Weyl candidate C a/ 3 y s- In our particular context we prefer to 
call this tensor the Bel- Robinson tensor of C a ^ in analogy with the Bel- Robinson tensor 
constructed outof the Weyl tensor [3]. The Bel- Robinson tensor of a Weyl candidate has 
the following properties — see [25j for detailed proofs of points (EJ), (EI]) and (En]); point 
( ITvj) follows by taking the covariant divergence of (13.81) and then using repeatedly (13.21) 
and fl3T3D . 

Theorem 2. IfT a pp U is the Bel-Robinson tensor of the Weyl candidate C a pp V then 

(y T{a^iiv) = T a /3p U , T alMJ = 0. 

(ii) Generalized dominant property: ifu^, u 1 ^, u^, are causal future- directed vectors 
then T a pp jV u1u2U l ^u v i > 0. This property admits an alternative formulation: let 
E = {eg, e^, e^, eg} be any orthonormal frame with the timelike vector. Then the 
following inequality holds 

^0000 > l^abcdl, 

where the component indices refer to the frame E. 
(Hi) T a pp V = <^=^> Cappy = <^=^> 3 a timelike vector u p such that T a/ 3p U u a u l3 u tJ, u 1 ' = 0. 
(iv) The covariant divergence of T a pp U is given by 

VqT^j = 4 C/3) a Scr — 2J^ aa C/3)Sacr ~ 9 'p-yJ^^ ^ 'saap- (3-9) 

Next, we explain how the properties summarized in theorem [2] are useful in the 
study of the causal propagation of Cp Vp \. We follow closely [4] in our exposition. Let 
D + {B) be the future Cauchy development of a closed achronal hypersurface B and pick 
up any event q G D + (B). Furthermore, let J~{q) denote the causal past of the event 
q. The set /C = J~(q) PI D + {B) is compact and it contains the set B = /C PI B. The 
boundary of /C is <9/C = B U H + {B) where H + {B) is the future Cauchy horizon of B. 
Also, since the interior of D + {B) is globally hyperbolic we deduce that it can be foliated 
by a family of spacelike hypersurfaces t G /, where to < t < t\, to < 0, < t± and 

Sq = B. The figure 13.11 shows a schematic view of this geometric construction. 
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Let n p be the timelike normal 1-form to the leaves of {Y, t }tei and define the quantity 

W(t)= [ T lu , afi 'nPn v n a n f) dM= [ dt' I [ T^ a ^n u n a n p dzA . (3.10) 

Jj-(E t )r\K Jo \jT, t ,n!C J 

Here, dM. is the positive measure on the spacetime M. constructed from the volume form 
V^pcr ( an Y measure defined from a volume form will be understood as positive) and d£ t ' 
is the measure on the hypersurface E t / obtained from its volume form n^^po-ls ,. The 
scalar T^ llJp xn fJ "n u n p n x is called the superenergy density and is everywhere non-negative 
— see point ([n]) of theorem [2] — which entails W(t) > 0. Also, from the definition of 
W(t) we get 

W'{t)= [ T^ aP nWn a n dZ t = [ P^ dZ t , (3.11) 

where P p is defined by 

P^ = T^ pa n v n p n' 7 . (3.12) 

Again point (Jn|) of theorem [2] implies that P^n p is also a non-negative quantity and 
hence W'(t) > 0. Now, according to the Gauss theorem we have 



dM= [ PX tffit - f PX dB+ [ P^ d(H + (B)), (3.13) 

JTittlK JB JH+(B) 



J-(£ t )n/C </£ t n/C JB JH+{B) 

where dB and d(H + {B)) denote the measures on the hypersurfaces B and H + {B) 
respectively induced by the volume forms n^^^lei an d ^V^pxIh+iB)- m ^ ms ^ as ^ 
expression k p is a causal vector which is orthogonal to the null hypersurface H + (B), 
points outward to /C (therefore it is future-directed) and is such that k fJ, r] IMUp x\ H +^ ^ 0. 
Note that these three conditions do not fix the vector k p univocally. By combining 
(13~TT|) and fl3TT3|) we deduce 

W'{t)= [ V p P p dM+ I P^dB- [ P p k^d{H + {B)). (3.14) 

Jj-(£ t )n)C JB ' JH+{B)r\J-(T, t ) 
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Now, we need to estimate the right hand side of (I3.14j) . According to point ([n]) of 
theorem [2], the quantity P p k^ is positive and therefore 



< / P% d(H + (B)). (3.15) 

JH+(B)nJ-(T, t ) 

On the other hand, if C a p^ v vanishes on B then W(0) = and 

P p n p dB = 0. (3.16) 



Next, we set up an orthonormal frame E = {eg, e±, e%, e^} with = and we consider 
the set of components V a n b in such a frame. Since these are continuous functions we 
deduce that in the compact set /C the estimate V a n b < mo, Va, b — 0, 1, 2, 3 holds for 
some constant rriQ. Hence 

V a P a = 3n c n d T cdab V a n b + n a n b n c V d T dabc 

3 

< 3m V T 00cd + n a n b n c V d T dabc < 3m T abcd n a n b n c n d + n a n b n c V 'jr*** , 



c,d=0 



(3.17) 



where in the last step, point (JD} of theorem [2] was used. From here, we deduce 



V^P" dM 



j-(£ 4 )n/c 



< 



I (V^T^n^pTia dM + 3m [ T^ p(T n^n v n p n a dM. (3.18) 



'J-(S0n/c JJ-(Y,t)c\K, 
Using the above estimates in (I3.14p we get 



< W'(t) < [ {V p T pvp °)n u n p n a dM + 



3m W(t). (3.19) 



Now it only remains to estimate n u n p n a V ^T^f" 7 and to that end (13.91) is used. The 
actual estimate will depend on the particular expression of the Weyl current. For 
example, the simplest case arises when J pup = where no further estimates are needed 
and Gronwall's lemma entails 

W(t) = 0, te[0,h), (3.20) 

which in turn implies that the scalar T pva fin p n v n a n^ vanishes in /C. From the 
arbitrariness of q we deduce that the superenergy density is zero in D + (B) and by 
a similar argument as the above one can show that this is also the case for D~(B). 
Hence point 3 of theorem [2] implies that C a p pu = on D ± {B). Consequently the tensor 
C a /3^ u propagates causally whenever the Weyl current is zero. 

Remark 1. Point ( lull) of theorem [2] entails C^^b = if and only if 
(T a/3fiu n a n l3 n fl n u )\i3 = 0. Accordingly, the causal propagation of a Weyl candidate 
can also be formulated in terms of the vanishing of the superenergy density at the 
hyper surf ace B. 
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4. The Mars-Simon tensor 

We start this section by recalling some basic background about vacuum spacetimes 
possessing a Killing vector — see e.g. [18J for a complete account. Let (M.^^) 
be a vacuum solution of Einstein field equations and let £ M be a Killing vector, 
V^;, + = 0. This condition enables us to define a 2-form, the Killing form, 

by 

F pv = V[^„] = V^. (4.1) 
Elementary manipulations lead to 

V p F uX = W vXtip ^. (4.2) 
The self-dual Killing form T pv is then 

J> = i^ + ii^, (4.3) 

with F* v = {l/2)r]^ Xp F xp the Hodge dual of F pu . From the previous definition it is 
straightforward to obtain the property (self-duality) 

= i •?>• (4-4) 
We shall write T 2 = T ilv T i> ' v . The Killing form is said to be degenerate if T 2 = 0. 
Other basic algebraic properties of T pv are 

V/ = ^^, V^ = ' V^ = 0- (4-5) 

From a Killing vector £ M we may define its associated Ernst 1-form by 

a u = 2^F pu = -V„A - ifl„, (4.6) 

where A = is the norm of the Killing vector and Q p = rj pi/Xp ^ v V x ^ p denotes its twist 
form. It is well-known that in a vacuum spacetime is closed — i.e. V^cr^j = 0. Thus 
locally — if we are in a simply connected region — there exists a scalar field a, the Ernst 
potential, such that o~ p = V^c The self-dual Weyl tensor and the self-dual Killing form 
are related by 

VpF vp = 2W up , a C, (4.7) 
a property which is easily obtained from (14.21) . Finally, we introduce the tensor X puXp 
defined by 

F pvXp = ^ {9p\9u P - 9p P 9v\ + vjiwXp) • (4.8) 

The tensor T puXp has the obvious symmetries T pvXp = F[p U ]x P = F Xppu . This tensor can 
be regarded as a metric in the space of self-dual 2-forms. 

We have now all the ingredients to define the Mars-Simon tensor. This is a four- 
rank tensor S pvXp whose expression is given by [jj 

<Spu\p = 2W^Ap + j — — i^F avFxp — QF 2 I puXp S j . (4.9) 

§ We follow the convention set in [19] of calling this tensor the Mars-Simon tensor. Note however, that 
our conventions are slightly different. 
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From the previous expression we deduce that in a vacuum spacetime admitting a Killing 
vector, the Mars-Simon tensor is determined up to a choice for the Ernst potential. We 
make a choice such that a ^ 1 on Ai. Our interest in the Mars-Simon tensor lies in the 
following result — see [2TJ [22] for a proof. 

Theorem 3. Let {M. ) g liu ) be a smooth, vacuum spacetime with a Killing vector I; 11 . 
Let T^v denote the associated self-dual Killing form. If there is a non-vanishing real 
constant M such that the conditions 

JF 2 = -M 2 (l - af, (4.10a) 
S^xp = 0, (4.106) 
hold on a non-empty M C M. then {J\f, g^v) is locally isometric to the Kerr spacetime. 
Condition (14. 10 aft of theorem [3] can be reformulated in the form 

S(M,O = 0, (4.11) 

where the scalar S(M, £) is defined for any Killing vector £ and any real constant M by 
the expression 

E(M,i) = T 2 + M 2 {l-af. (4.12) 
A very important property of the scalar H(M, £) is that it fulfils the differential equation 

= -2T^S waP e + (4.13) 

To see it, differentiate both sides of (14. 12p and use (14. 7p to replace the covariant 
derivatives of T^ v . The result is 

= -4(J^W MpajB e p + M 2 a,(l - af). (4.14) 

Next use ( 14.91) to write W MpQ /3 in terms of S wa p and replace all the occurrences of 2^ pal3 
by means of (14.81) to obtain 

V M 5 = 2(2MV> - I) 3 - ^S, paP i") + 1 » P \ ^ -(4.15) 
We use here ( 14.41) and (14. 6 p with the result 

V^-2^^+ 4 ^ 2 + A ?- j)4) . (4.16) 

a — 1 

from which ( I4.13P is deduced. 

4-1. Mars- Simon tensor as a Weyl candidate 

In [19] it is proved that the Mars-Simon tensor S^ap is a Weyl candidate. Furthermore, 
it was shown that its Weyl current is a linear function of S^ap in which case the wave 
equation (13. 4p is homogeneous in S pua p. The derivation of these results are discussed as 
some of the intermediate calculations will be needed in section [5j 

Proposition 1. The Mars-Simon tensor is a Weyl candidate. 
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Proof. From (14.91) it is clear that S^ aX = <S^[ a \] = <Sa\pv Using (14. 5 p one can show 
that S „ p p \ = and therefore the Mars-Simon tensor is traceless. To finish our proof 
we need to show that the Mars-Simon tensor fulfils the cyclic identity. To that end we 
compute S^px] getting 

6 i T 2 

Sp[v P \] = YZT^^^p^ ~~ 2(1 - a )^ pA ' ( 4 ' 17 ' ) 

The right hand side of this expression vanishes by virtue of the identity 

FulvFpX] = y^i F 2 Vp»p\, (4-18) 

which in turn is a consequence of the algebraic property Tp^Tpx] = F\ pv T p x\ and the 
fact that in four dimensions J-\p, v 3~p\] is proportional to r] pup x- □ 

A direct computation shows that the Mars-Simon tensor is self-dual 

<W = i S; vpX . (4.19) 
Therefore Re(S pup \) and Iul^^a) are then related by duality as explained in section 

m 

lm(S pupX ) = ^(5/), Re(S pupX ) = ~ VpXa plm(S^). (4.20) 

Any property of S pup \ will admit an equivalent counterpart formulated in terms of either 
its real or imaginary part. In particular, this shows that both ~Re(Sp, up x) and Im^S^px) 
are Weyl candidates and any of them contains the same information as the original 
Mars-Simon tensor. 

Proposition 2. The covariant divergence of the Mars-Simon tensor is given by 

V a <5 a w = — — - {^{T^S^pap + T^SySap) + Zg^SsiaXpF^ C a - (4-21) 

Proof. To calculate V a S a i 3lS we start from (14.91) and take the covariant divergence of 
S pu xp- Next we apply the relations (g2D, V^cr = a p , V a W a apu = and V a l pupa = 0. 
Finally we expand all the occurrences of X pvpa using ( 14. 8 \ with the result 

V a S a ^ s = - —^ {-YlFpuT^a 01 + (i f]f3js a o- a + 2g f3 y 1 a S \)J r2 ) 

+ ^((21 m,s a w; nK e + Agp [s w 1]a u a )F™ - i2^/w 7foP r) (4.22) 

The dependence on the Ernst 1-form can be removed from this expression if we write 
a p in terms of T pv and ^ by means of (14.61) . The final step is to rearrange the terms in 
the resulting expression in order to render it in the form of (I4.2ip . This is accomplished 
by using (14.91) to write V\? pv x P in terms of the Mars-Simon tensor, expanding again all 
the occurrences oiX pvpa by means of (14.81) and using (14.51) . (14.41) where necessary. □ 
Equation (14.211) enables us to define the Weyl current of the Mars-Simon tensor as 



J(i-y& = — 

a 



— j" {^{F[5 P S-y\flap + F/S-fSap) + f3[iS 5}aa pF* P ) £° (4.23) 
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Since j7/? 7< s is linear in the Mars-Simon tensor, we conclude — in view of (I3.4p — 
that S pvpa fulfils a homogeneous wave equation. This interesting property was used 
in [19] to show that the domain of outer communication of an asymptotically flat 
smooth stationary spacetime is, under certain additional conditions, isometric to the 
Kerr domain of outer communication — i.e. the uniqueness of Kerr. 

Using V a S a pup = J pvp together with the self-duality of the Mars-Simon tensor one 
deduces 

3pvp = i <J pV p = ~^}vpap3 p ^ i (4.24) 
that is to say, J pvp is also self-dual with respect to the block of indices up. 

4-2. Orthogonal splitting of the Mars- Simon tensor 

In our forthcoming calculations it is necessary to use the orthogonal splitting of S pup \ 
with respect to a unit timelike vector field n M . Since S pvp \ is a Weyl candidate we can 
use (12.121) to find such a orthogonal splitting. The resulting split is given by 

Sp, vp \ = 2(^ A 7^ M + lp[ p T\]v + i £p\ a n[pT u \ ° + i ^pva^pT^ ")> (4-25) 

where 

T pu = S ppvX n"n x . (4.26) 

The tensor T pu is symmetric, spatial and traceless. Alternatively, the orthogonal 
splitting of the Mars-Simon tensor can be calculated directly from (14.91) . To that end 
one needs to find the orthogonal splitting of the different quantities which appear in 
(14. 9p . The orthogonal splitting of W pup \ is given in (12 . 1 2[) and the orthogonal splitting 
of T pv is calculated by means of (12. 8j) and is given by 

^> = i e ppv 8 p - 2£ [lx n u] , (4.27) 

where 

£ p = T p y. (4.28) 

One also needs the orthogonal splitting of I pvp \ which is easily calculated from (14.81) 
and (12.51) so that 

Ipupx = h^phxju + 2n[ u h p ][ X n p ] + i (n[ p e u ] pX + ri[ p ex\^)- (4.29) 

We insert relations (I2TT2]) . fT4~2TD and ffl~2"9l into the formula (EDI and then equal the 
resulting expression to (14.251) . The result of that is the relation 

T pv = 2£ pu + (3£ p £ u - h pu £ a £ Q ), (4.30) 

1 — 0" 

which shall be used later on. Also needed later on is the relation 

T 2 = -A£^ =^ S(M, = -4^£" + M 2 (l - a)\ (4.31) 
which is derived from equation (I4.27p . 
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The Killing vector ^ can be decomposed in the form = — Yn p + Y p where Y, 
Y^ are called the Killing lapse and the Killing shift respectively. We need to find a 
formula relating £ p to Y, Y p . To that end, we start by taking the covariant derivative 
of both sides of £ M = —Yn p + Y p and replace V p n u by means of (12. 4p . The result is 

= Y(K^ V + A v n^) + VJfy - n u X7 p Y. 

Next we use in this expression the property V p Y = D p Y — n p L^Y and equation (12.141) 
to replace V P Y and V P Y V respectively yielding 

= YK^ + D p Y u + n p n u (C H Y - A"Y p ) 

+ n„(YA u - K vp Y p - C H Y U ) - n u (K w Y» + D^Y). (4.32) 

The antisymmetric part of this expression gives the orthogonal splitting of the 2-form 
which is 

= -YA [lx n v] + n [p D u] Y + D [p Y u] - n^C^y (4.33) 

Also the symmetric part of (I4.32p is the Killing condition V( At ^ i/ ) = 0. Explicitly one 
has 

YK^ + D {p Y v) + n^CnY - A P Y P ) + n (lx (A u) Y - 2K u)p Y p - D v) Y - C H Y u) ) = 0, 
from which we deduce 

£ H Y p = A^Y - D,Y - 2K lxp YP, (4.34) 
which renders the relation (14.331) in the form 

Ffj, u = 2n ]/l K v]t> Y" + D [IM Y U] + 2n [fM D u] Y. (4.35) 
From this expression we may compute the orthogonal splitting of F* v which is 

F; u = e^KjY" + D?Y) + n [p e u]a(3 D a Y^ (4.36) 

where equation ( 12. 5ft was used along the way. Inserting equations ( 14.351) and ( 14.36[) in 
( 14. 3 p and applying ( 14.281) on the resulting expression we deduce 

= K, p Y p + D p Y - h e^D'Y", (4.37) 
which is the required relation. 



5. The causal propagation of the Mars-Simon tensor 

In this section a proof of the following result is provided. 
Theorem 4. The Mars-Simon tensor S pup \ propagates causally. 

Proof. According to ( 14.201) it is enough to show that either Re(S pup \) or Im^S^px) 
propagates causally. We choose to work with the former, so let us set S pvp \ = Re(iS M „ P A) 
and define L^x = Re(J' IMU x). Equation (14.211) entails 

V a S a ^ s = Lp lS . (5.1) 
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Hence, Lp^s is the Weyl current of the Weyl candidate S^px- To study the causal 
propagation of S^ vp \ we follow the general procedure explained in section 13.11 We 
denote by B pvp x the Bel- Robinson tensor constructed from the Weyl candidate Sp Vp x- 
To prove our result, we need to find a good estimate for the quantity n u n p n x \7 pB^ vpX 
— see equation (13. 19ft . One can write 

IT^S^x + 4L»;n?n x n a Sx, a v, (5.2) 

where ( 13.91) was used. In order to work out the right hand side of equation (I5.2p we 
calculate the orthogonal splitting of S^pX and L^x with respect to n M . The tensors 
Spupx and L pup are related to Sp Up x and J pvp by means of the relations 



n»n»n x V p B» upX 



S. 



fiupcr 



2 {S pupa 



+ S 



livpa ) i 



pv p) i 



(5.3) 



and therefore their orthogonal splittings can be calculated once those of S pvpa and 
J pvp are known. The orthogonal splitting of S pvpa is given in equation (I4.25P and the 
orthogonal splitting of Jp Vp can be calculated from equation ( 14.231) given that we know 
the orthogonal splitting of all the quantities which appear in the definition of Jp, vp — see 
(lA.lj) and (1A.2j) - (lA.5l) in appendix Appendix A.l for the precise formulae. Inserting 
these orthogonal splittings in (15.21) we get 



P n X VpB» upX 



12 Re 



a 



- 6i e^T^T/Y 



£' 



a 



1 



a 



1 



(5.4) 



Now, using the equations (1A.13I) - (1A.14I) of appendix Appendix A. 2 we can deduce the 
property 

B a p U p(£ a n»n p Y^£ a n»n v nPY) = -i e^T^T p a Y+£»Y,T^T aP -2£ »Y v T {v 



which when combined with ( 15. 4p yields 



n»n^V p B» upa 



Re 



12 



—B pua g(£ tM n f3 n u n a Y - £»n a n Y l 
a 



(5.5) 



Equation (15.51) is the key to get the required estimate. To obtain it, let us introduce an 
orthonormal frame E = {e^, e^, , } with = . Since £ M , Y^ are spatial vectors 
they can be written as a linear combination of the vectors , , eg . Therefore, point 
(Jul) of theorem [2] entails the estimates 

12 



Re 



B pvoiP £^n»n a Y 



-Re 



- a 
12 



a 



-B^ a ^n a n p Y v 



< m l B llua!3 n IJ 'n u n a n 13 



< m 2 B puaP n^n y n a n 13 



(5.6a) 



(5.66) 



for some scalar continuous functions m\, m 2 defined on JC. Using these estimates in 
equation (15. 5p we obtain our final estimate n v n p n a \7 pB^ vpa < m^B PiUa pn^n v n a n 13 where 
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m 3 is the upper bound of m x + m 2 in the compact set /C and thus (13.191) can be written 
in the form 

< W'{t) < (m 3 + 3m )W{t), Vte[0,ti). (5.7) 

Since W(0) = 0, then an application of Gronwall's lemma enables us to conclude that 
W{t) = 0, Vt G [0,ti) which, according to the considerations made at the end of section 
13. II proves the causal propagation. □ 

Remark 2. We stress that the fact that S pup \ satisfies a wave equation and the property 
(I4.2ip do not, by themselves, suffice to show the causal propagation of this tensor. If we 
impose that S pup \ and some directional derivative of this quantity vanish on a spacelike 
hypersurface B then using the wave equation and standard results on hyperbolic partial 
differential equations allow us to conclude that S pvp \ vanishes on a neighbourhood of B. 
In doing so, extra conditions needed to be imposed so that S pvp \ is indeed zero in D{B). 
This approach leads, in principle, to a weaker result than the one given in theorem HI 

Consider now any vacuum spacetime admitting a Killing vector £ M and use £ M to 
construct the Mars-Simon tensor S pvp \ and as above, define its real part S pup \ and the 
Bel- Robinson tensor B pvp \ of S pup \. For any unit timelike vector n M , define the positive 
scalar 



B pvpa n^n p n a . 



$(n,0 

One can find an explicit expression for £) if we use the property 
f TAl2l) in 



Appendix A. 2 



B 



[ivpa 



rfn v n p n a 



T T 

1 pv £ 



■pv 



4£^£^-24Re 





-12 




2 

+36 


£»£» 




a - 1 




(7-1 






a - 1 



(5.8) 

-see equation 



(5.9) 



which shows that $(n, £) can be expressed in terms of S pu , S p and a exclusively. When 
$(n, ^) vanishes on a Cauchy hypersurface, the causal propagation of S pvp \ entails the 
following result. 

Theorem 5. Let {M. 1 g pi ) be a globally hyperbolic vacuum spacetime and suppose that 
M. admits a Killing vector If there are a Cauchy hypersurface B C M. and a real 
constant M such that $(iV, = and 5(M, = where N is the unit timelike 
normal to B, then M. is locally isometric to the Kerr spacetime. 

Proof. If $(iV,|) is zero on B then by ([EH]) B,j, vpa N^N v NPN a \ B = 0. Therefore, 
remark [1] leads us to 

S p upa \b 0. 

Using that the Mars-Simon tensor propagates causally, one deduces that S pupa = on 
M.. In particular, this implies that (I4.13P is rendered in the form 

(7—1 

which can be integrated to yield S = A(l 
condition H(M, = implies that A 



(5.10) 

a) 4 for some complex constant A. The 
and hence H (M, ^* ) vanishes on .M . 
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Theorem [3] can now be applied and we conclude that M. is locally isometric to the 
Kerr spacetime. □ 

6. Application: characterization of Kerr initial data 

We explain in this section how to use previous results to find conditions which ensure 
that the data development of a vacuum initial data set is isometric to a subset of the 
Kerr spacetime. The basic idea behind is to express the result of theorem [5] in terms of 
conditions on a vacuum initial data set (£, hij, Kij). 

Theorem 6 (Kerr initial data). Let (£, h^, K^) be a vacuum initial data set and 
assume that there exist two scalar fields Y , a, a vector field Yj and a real constant M, 
all defined on S ; fulfilling the following conditions 

£j£ j = ]m 2 {1-&)\ (6.1a) 

2 



4 



AS^Sij - 24Re 



a - 1 



- -M 4 |l -a| 6 + 36 
4 1 1 



£ % £i 



GT- 1 



(6.16) 



Dja = 2YK jl Y l - 2 i s jm i(K k l Y k Y m + Y m D l Y) + TYDjY + 2Y l D {l Y jh (6.1c) 
D(jYj) + YKij = 0, (6.1 d) 

DiDjY + CytKtj = Y{r i3 + KKy - 2K iX K\), (6.1 e) 

where 

Ej = K jk Y k + DjY - h e ]kl D k Y\ (6.2a) 

E jk = ^(E jk -iB jk ). (6.2b) 

Then the data development (M, g^ v ) hij, K^) is locally isometric to an open subset 

of the Kerr spacetime. 

Remark 3. The conditions (16. Id)) and (16. lei) are the Killing initial data (KID) 
equations. It is well-known that if both a solution (Y, Yj) and the initial data set 
(S, h^, Kij) are suitably smooth, then the development of the initial data (E, h^, Kij) 
possesses a Killing vector such that the pull-back of its restriction to S coincides with 
the given Killing data — see for example, P, [13] , and in particular [8] for an explicit 
set of smoothness conditions on the Killing data. Alternatively, if the Killing data is 
transversal, that is Y ^ 0, then one can make use of the concept of Killing development 
to ensure the existence of a spacetime with a Killing vector associated to the Killing data 
under question — see p]. Again, under suitable regularity conditions it can be shown 
that the Killing development always includes the Cauchy development of (£, h^, K^) 
— see again (9): the Killing development is unique and maximal among the class of 
spacetimes containing the relevant Cauchy data and a Killing vector with non-closed 
orbits; hence, if the Cauchy development has a timelike Killing vector, then it must be 
a subset of the Killing development. 
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Proof. Let M. be the data development of the vacuum initial data set (E, hij, Kij) and 
let £ M be the Killing vector mentioned in remark [3) As usual consider a foliation 
t G /, of M with So = and let n p be the unit timelike normal to the leaves of 

this foliation. Let us adopt the formalism and notation explained in section H] with £ M 
being the Killing vector. Then the orthogonal splitting of is = —Yn p + As 
explained in remark [3] we have the property 

<f)*Y = Y, <9^% = Yj. (6.3) 

which combined with equations 1 16. 2 qj) and (14.371) entails dj(p^(£^) = £j. The orthogonal 
splitting of is calculated on one hand from equation (14.61) using the expressions (I4.27p . 
(T437]) and the property ^ = -Yn^ + Y*. The result is 

0> = 3^ + sn^, (6.4) 

where 

s = -2K^Y V - 2Y"D tt Y + i e wv Y»D p Y v , 

s M = 2YK W Y" - 2 i e flup (K a p Y a Y u + Y v D P Y) + TYD^Y + 2Y v D [v Y lA . (6.5) 

On the other hand, the orthogonal splitting of <r M can be calculated from the relation 
&n = V M a yielding 

s = -Cfta, Sf, = D^a. (6.6) 

Combining the previous equation with the pull-back of equation ( 16. 5ft under <fi and (16. left 
we deduce 

Dj{a) = DiWip)). (6.7) 

Upon a suitable choice of the Ernst potential a, this equation can be integrated on E 
to give </>*(<t) = o. Next consider the Mars-Simon tensor S^ U p\ constructed with the 
Killing vector and calculate its orthogonal splitting with respect to n p — equation 
(I4.25p . Using 114. 30p and our previous results, we get 



(6.8) 



Expanding this product and using conditions fl6.lap - 116.lfrl) we conclude 

*(n,|5U (E) = 0. (6.9) 

Moreover, equation (16. laf) implies H(M, OU(E) — — see (14.3 ip . Thus, theorem [5] applies 
and therefore we conclude that (A4,gp U ) must be locally isometric to an open subset of 
the Kerr spacetime. □ 
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6.1. Schwarzs child initial data 

Theorem [6] asserts that under certain conditions the development of a vacuum initial 
data set is isometric to an open subset of the Kerr spacetime but nothing is said as to 
the resulting development being one of the specializations of the Kerr spacetime: the 
Schwarzschild solution. In this section we show the form which theorem [6] takes in that 
case. This complements previous work of the present authors about the same issue [16J. 
To state our result we need a preliminary lemma whose proof can be found in [6]. 

Lemma 1. Let (E, hij, Ky) be a vacuum initial data set and assume that there exist a 
scalar Y and a vector field Yj, all defined on E ; fulfilling the following conditions 



Then there exists an integrable Killing vector £ M in the data development of (E, hij, Kij). 

With the aid of this lemma we can now derive conditions guaranteeing that 
the development of a vacuum initial data set is isometric to an open subset of the 
Schwarzschild spacetime 

Theorem 7 (Schwarzschild initial data). Let (E, hy, K^) be a vacuum initial data 
set and assume that there exist two scalar fields Y , b and a vector field Yj, all defined on 
E, fulfilling the conditions of theorem^ with Ii6.1d\) - /[6.1e\) replaced by Ii6.10a\) - ll6.10b\) . 
Then the data development M. is locally isometric to an open subset of the Schwarzschild 
spacetime. 

Proof. The first step in the proof of this theorem is to use lemma [1] to show that the data 
development Ai admits an integrable Killing vector ^ . The remaining part of the proof 
follows the same pattern as that of theorem [6] and therefore we reach the conclusion that 
M. is isometric to an open subset of a specialization of the Kerr spacetime possessing an 
integrable Killing vector which, as is well known, is the Schwarzschild spacetime. □ 

Remark 4. It is worth recalling that the analysis of [16] has provided an explicit formula 
for a timelike KID candidate in terms of concomitants of the initial data (hij,Kij). 
Hence, a combination of theorem [7| with the such timelike KID candidate renders an 
algorithmic characterization of Schwarzschild initial data alternative to the one given in 



7. Conclusions 

The main result obtained in this article, theorem [6] provides a characterization of 
initial data for the Kerr spacetime by first requiring that a set of overdetermined 
elliptic differential equations — the KID equations — admit a solution; and then in turn 
requiring that certain objects constructed out of the solutions to the KID equation 
satisfy very particular relations. This characterization fails to be algorithmic in two 
ways. First, it is in practise a non-trivial problem to verify whether a certain initial 



YD[iYj] + 2Y [i D j] Y + 2Y [i K ii iY l = 0, 
Y {l DjY k] = 0. 



(6.10a) 
(6.106) 



1161. 
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data set admits a KID or not. A pointwise procedure for doing this has been discussed 
in [2]. Alternatively, one could resort to a non-time symmetric generalization of the 
geometric invariant characterizing static time symmetric initial data sets constructed 
in [10] — which requires solving a fourth-order elliptic partial differential equation. In 
any case — and this leads to the second way the result fails to be algorithmic — even if 
one could ascertain the existence of a KID on the initial data, one requires to know the 
KID in an explicit manner in order to be able to verify the conditions (I6.1al) . (I6.1&I) and 
( 16. left . It is plausible that from an hypothetical characterization of the Kerr spacetime in 
terms of, say, concomitants of the Weyl tensor — analogous to that for the Schwarzschild 
spacetime given in [1?>] — one would be able to construct, following the ideas of [16], KID 
candidates for the initial data set (£, hij, Kij). If such a KID candidate were available, 
then theorem [6] would constitute a powerful tool to analyse the behaviour of dynamical 
black hole spacetimes, both from an analytic and a numerical point of view. 
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Appendix A. Technical details about the calculations 

Appendix A.l. The orthogonal splitting of Jp Vp 

Given the symmetries of the Weyl current Jp Vp it is immediate to realise that its 
orthogonal splitting must take the form 

Jiwp = ripZ up + n [p X u] p + npn [u Vp] + Vp Vp . (A.l) 

It takes a greater effort, however, to find the explicit expression of the spatial tensors 
Zp V) Xpv, Vp, Vp Vp . These expressions are found by inserting into equation (14.231) the 
orthogonal splitting of Sp Up \ shown in (14.251) . the orthogonal splitting of Tp V shown in 
(14. 271) and the relation = —Yn p + Y 7 ' . A computer calculation yields the results 

Zp V = ^—(i YZ a \ p e u]ap + 3£ { p% ]p Y p + Y [u Tp ]a £ a - i Ep Vp £ a T a P Y), (A.2) 
X = - ( £ a (AiY(X p f \ - f T p ) + 6T Y - 2T Y - AT Y + Ah T Y p ) 

< v pv — -y y-' ^ J V (ot u )l- i P c vap- 1 p J ~ u ■'■up 1 a ^ - 1 pa 1 u ^■'ua 1 p ~ ^"'up * ap 1 ) 

+ 6£ u Tp a Y a )y (A.3) 
Vp = ^ (4i Y p {Sp Kp T a) A + e apX T^ A ) - 2YT„ a ) , (A.4) 
Vp Vp ee ^— \2 i e upX £pT a x Y a + 2£ [p (3YT u] p + i e v] p X T a x Y a ) 
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+ £ a (^ — 2\(e paX T\ u x F M ] + e uaX T[ p x Y p } + 2e paX T[ u x Y p ] 

Y [B upa T pX -\- Sfia^i/T 'p]; 



(A.5) 



Appendix A. 2. The orthogonal splitting of B pupa 

In [Hj the orthogonal splitting of the Bel-Robinson tensor was calculated and the 
different parts resulting from such splitting were studied in detail. The calculations 
are similar for the Bel-Robinson tensor of a Weyl candidate and therefore we only need 
to adapt these results for the particular case of S pvpX . In this case the orthogonal 
splitting of B pvpX reads — see section (5.1) of [14] — 

B pvpX = Wn p n u n p n x + AP {p n u n p n x) + 6t^ u n p n X ) + 4:Q iflup n X ) + t pupX , (A. 6) 

where the different parts of this splitting are given by 

W = e a pe aP + b a pb^, (A.7) 
P u = 2b a x e pX e p ap , (A.8) 
V = Wh pu - 2{bJ>b vp + e/e up ), (A.9) 
Qvwp = h pu P p - 2 {b ua e pP + b pa e ul3 ) e p al3 , (A. 10) 

tpvp\ = ^( b nv b P \ + e pu e P \) ~ hp\tfiv + 2h(\tp)p + 2h^xt p ) v - h pu t pX + 

+ W{h pv h pX - 2h p[p h x)v ), (A.ll) 

where e pu and b pu are respectively the electric and magnetic parts of the Weyl candidate 
S puXp . From (tOI) and (fA~T0l) we get 

B pupX nVn p n x = T pX T pX (A. 12) 

B pupX n»n?n x = -n p T px T pX + i e pvp T Xv T x P . (A. 13) 

B pvpX n p n x = -2\ p T u)p + (n p n u + h pu )T px T pX - 2i T ap T a X e pX(u n u) , (A. 14) 

where the relations 

&p,u Tpu ~\~ 1 p,wi b pv i {T~ P u T P u) i 

need to be used. 
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